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Abstract 


Summary 


Max, Timothy A.; Souter, Ray A.; O’Halloran, Kathleen A. 1990. Statistical 
estimators for monitoring spotted owls in Oregon and Washington in 1987. Res. 
Pap. PNW-RP-420. Portland, OR: U.S. Department of Agriculture, Forest Service, 
Pacific Northwest Research Station. 13 p. 


Spotted owls (Strix occidentalis) were monitored on 11 National Forests in the Pacific 
Northwest Region of the USDA Forest Service between March and August of 1987. 
The basic intent of monitoring was to provide estimates of occupancy and reproduc- 
tion rates for pairs of spotted owls. This paper documents the technical details of the 
statistical estimation procedures used to produce all the estimates resulting from that 
effort. All relevant formulae are presented as well as decisions made about collaps- 
ing strata to estimate variances. Some complicating aspects of the monitoring effort 
are also discussed. 


Keywords: Spotted owl, Strix occidentalis, sampling, stratified random sampling, 
monitoring. 


Spotted owls (Strix occidentalis) were monitored on 11 National Forests in the Pacific 
Northwest Region of the USDA Forest Service between March and August of 1987. 
The objective was to provide estimates of occupancy and reproduction rates for three 
distinct subpopulations consisting of different land classification types as defined in 
the Pacific Northwest Region. Over time, point estimates from 1987 will be used with 
point estimates from other years to establish trends. The sampling design consisted 
of stratified random sampling independently designed and implemented for each 
subpopulation. This paper documents the technical details of the statistical estimation 
procedures used to summarize the data obtained from sampling. All relevant formu- 
lae are presented as well as decisions made about collapsing strata to estimate 
variances. Field protocols used to obtain all measurements may have resulted in the 
biasing of some estimates. Two sources of potential bias were misclassification of 
occupancy of sites and failure to detect multiple occupancy. These sources of poten- 
tial bias are thoroughly discussed. 
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Introduction 


Population 


During summer 1987, spotted owls (Strix occidentalis) were monitored on National 
Forest land in the Pacific Northwest Region (Region 6) of the USDA Forest Service. 
The general objective of monitoring was to provide estimates, for three distinct land 
classifications, of occupancy by pairs and of reproductive rates of the spotted owls. 
The results of this monitoring as well as details about the objectives and methods 
were provided in a report prepared by the USDA Forest Service.’ The purpose of 
this paper is to document the technical details of the statistical estimation procedures 
used to produce the results presented in that report. 


To understand the estimation procedures used, it is necessary to have a clear under- 
standing of the three distinct land classifications used in the sampling process. These 
three land classifications, defined completely in the next paragraph, specify three dis- 
tinct subpopulations of land area within Region 6. Each of these subpopulations was 
sampled independently by using area sampling methods. It is important to clearly 
understand that in this sampling design, the term “subpopulation” refers to a distinct 
type of land area defined by the land classification definitions given in the next 
paragraph. 


Definitions of land classes were based on two factors. First, land was classified as 
“reserved” or “nonreserved.” Reserved land is specifically designated as “withdrawn 
from timber harvesting activities,” as specified by law or administrative decision by 
the Chief of the Forest Service; examples are wilderness areas and research natura! 
areas. Nonreserved land is generally available for timber harvesting. This class 
included some land currently classified as “not suitable for timber harvest.” The 
second factor was whether or not land was designated as a spotted owl habitat area 
(SOHA).7 In 1987, SOHAs were established on both reserved and nonreserved land. 
Thus four distinct land classifications were reserved non-SOHA, reserved SOHA, 
nonreserved non-SOHA, and nonreserved SOHA. Reserved non-SOHA land was not 
specifically included in the population for this monitoring effort. The three subpopula- 
tions sampled in 1987 were, therefore, the other three land classifications; that is, 
reserved SOHA, nonreserved non-SOHA, and nonreserved SOHA. 


Within the three defined subpopulations, land area was further classified into physio- 
graphic provinces (Franklin and Dyrness 1973) and National Forests. Eleven National 
Forests were included in the 1987 monitoring, because as of March 1987, SOHA 
networks had been well established on 11 of the 13 National Forests in Region 6 
where spotted owls occur. The 11 National Forests and physiographic provinces are 
identified in table 1. 


‘USDA Forest Service, Pacific Northwest Region. FY87 an- 
nual report: spotted owl inventory and monitoring. Portland, 
OR. Administrative report. On file with: Pacific Northwest 
Region, P.O. Box 3623, Portland, OR 97208-3623. 


? Guidelines used for establishing and spacing SOHAs are 
contained in a letter from the Regional Forester to Forest 
Supervisors: letter 1920, February 9, 1983. On file with: Fish 
and Wildlife Staff, Pacific Northwest Region, P.O. Box 3623, 
Portland, OR 97208-3623. 


Table 1—Physiographic provinces, acronyms used to reference the 
provinces, and National Forests monitored for spotted owls in 1987 


Province 
Physiographic province acronym National Forests 
Olympic Peninsula OLPE Olympic 
Washington Cascade Range WACA Mount Baker-Snoqualmie 
Gifford Pinchot 
Oregon Cascade Range ORCA Mount Hood 
Willamette 
Umpqua 
Rogue River* 
Deschutes 
Winema 
Coast Ranges of Oregon CORA Siuslaw 
Klamath KLAM Siskiyou 
Rogue River? 


2 The Rogue River National Forest was divided between the ORCA and the KLAM provinces. 


Most estimates were provided separately for the three defined subpopulations. Some 
estimates were provided by physiographic province, and others were made for groups 
of physiographic provinces. Other estimates were provided only for the specific sub- 
population for the entire portion of Region 6 that was monitored. No estimates were 
provided for individual National Forests. 


Sampling Design The sampling design used in each of the three subpopulations was stratified random 
sampling. Strata were the National Forests. Sample sizes were independently estab- 
lished for each of the three subpopulations as described in the “FY87 Annual Report” 
(see footnote 1). Within subpopulations, sample units were allocated among the 
National Forest strata approximately proportionally to stratum size. 


Different sampling units were defined for the different subpopulations. For the 
reserved and nonreserved SOHAs, sampling units were the 1000-acre designated 
SOHAs. The SOHA definition used for 1987 monitoring was 1000 acres of suitable or 
capable spotted owl habitat delineated within a 1.5-mile-radius circle (4500 acres) 
centered on a location of owl activity. Forests had delineated 2200-acre SOHAs for 
planning purposes, but 1987 monitoring was based only on the field work conducted 
within the 1000-acre boundary. In the nonreserved non-SOHA land, sample units 
consisted of all the area in a 1.5-mile-radius circle, or 4500 acres. 


Estimators for 
Stratified Sampling 


Notation 


Case 1, Continuous 
Variables 


Three basic types of stratified estimators were used in this survey, depending on the 
type of response variable measured. The three types of response variables were (1) 
continuous variables, (2) simple categorical variables with two possible classes, and 
(3) more complicated categorical variables. This more complicated categorical vari- 
able occurred when an observation was first classified as a simple categorical vari- 
able and then one of the original classes was further subdivided into two possible 
categories. Appropriate estimators will be presented for each of these three types of 
response variables. 


The following notation will be used throughout. This notation and the formulae used 
follow Cochran (1977, chapter 5). 


his a subscript denoting a stratum and ranging from 1 to L, 
Lis the number of strata in a population (that is, the number of National Forests in 
the population), 
i is a subscript denoting a particular sampling unit within a stratum, 
Nh is the number of units in stratum h, 
LL 
N is the number of units in the population = > Nh, 
h=1 
Wh = Nb/N is the stratum weight, 
Nh is the number of units selected in the sample from stratum h, 
L 
nis the total sample size = > nh, and 
h=1 
fh = Nh/Nh is the sampling fraction in stratum h. 


If the response variable was continuous, then the straightforward stratified estimator 
for the mean is 


Yst = » Whyh, (1) 
h=1 


L 
where 
Nh 
Yh =(1/nh) %  yhi, and 
i=1 
Yhi = response measured on the ith sample unit in stratum h. 


The estimated variance of Yst is 


L 
VA) = OS WK GSR) sicin- 
h=1 


Case 2, Simple 
Categorical Variables 


1 Nh 
sh? = —— > (yhiyh)* is the sample variance for stratum h. 
Nh-1 i=1 


The standard error of yst.is 
e. (Yst) = [v(¥sn] . (2) 
Approximate 100(1—c«) percent confidence intervals for yst are constructed as 
stt(t)[S.e.(Yst)] , (3) 
where 
t = the o/2 percentage point of the t-distribution with ne degrees of freedom. 


L 
ae (> QhSh*] al Qh?sh4/(nh—-1)] , (4) 


where 
Qh = Nh(Nh — nh)/nh , 


and is an approximate method of assigning an effective number of degrees of 
freedom to s.e.(yst) (Cochran 1977, equation 5.16). The construction of confidence 
intervals by this method assumes Yst is normally distributed. The computation of ne 
assumes the variable, yni, is normally distributed. 


If the response variable is categorical, with two possible responses, then the objec- 
tive is to estimate the proportion in one of the classes. To develop an estimator for 
this situation, define an as the number of sample units in stratum h with the desired 
characteristic (that is, the number of sample units in the class of interest). 


LL 
Pst = > Whbh, (5) 
h=1 
where 
Aw 
Ph = ah/Mh . 


Case 3, Complicated 
Categorical Variables 


: . 6 LoS ep 
An unbiased estimator of variance of Pst is 


L 
vipsy = © Wh?[v(pn)] , 


h=1 
where 
V(Ph) = (1-fh) (Bh) (1-Ph)/(nh-1) . 
The standard error of Dst is then 
s.e.(Pst) = [v(pst)] 2 (6) 


Approximate 100(1-a) percent confidence intervals can be computed for Det as 
Pst + t[s.e.(Pst)] , (7) 

where 

t is the o/2 percentage point of the t-distribution with ne degrees of freedom, and 

Ne is computed as in equation (4) with 


Nh 

2 ES AW AW 

Sh — (Ph)(1-Ph) and 
Nh-1 


Qh is as defined before. 


Although these approximate confidence intervals can be computed, the necessary 
assumptions are not met. The necessary assumptions are that Pst is approximately 
normally distributed and, to compute Ne, that the variable, yni, is normally distributed. 
Clearly yni is a categorical variable with values of either 0 or 1 and does not meet the 
assumption of normality. So the nominal probability level of these confidence 
intervals will be smaller than the actual probability level. The amount of the difference 
between the nominal and actual probability levels is not known. 


For complicated categorical variables, sample units were initially classified into one of 
two classes as for case 2. We let Ao and Aj represent the two classes. After the 
sample was Classified by this criteria, units classified as Ai were examined further 
and classified by another criteria into one of two categories, say Bo or Bi. The 
objective was to estimate the proportion classified as Bi within those classified as 
Ai. An example of using this more complicated categorical variable was in estimating 
reproduction. We let the A criteria represent occupancy by a pair; that is, units 
Classified Ao were not occupied and units classified A1 were occupied. We let the B 
criteria represent occurrence of reproduction; that is, units classified Bo had no young 
owls observed whereas units classified B1 had young owls observed. One objective 
was to estimate the proportion of occupied sites where reproduction occurred. 


Estimating 
Occupancy 


First, we will define some notation to help understand this situation. Let an be the 
number of units classified Ai, and bh be the number of units in Aj that are classified 
Bi. 


Within a particular stratum, the estimated proportion of A1’s classified as Bi is bh/an. 
And a consistent estimate of the variance of bh/ah (from Sukhatme and Sukhatme 
1970: 27-30) is 


v(bh/ah) = (1—fh) (bn/an2) (1— bh/an) - 8) 


Putting these estimators into the stratified estimators results in an overall estimate, 
for the population, of the proportion classified as A; that are also B1. We will call this 
estimated proportion (p11)st. So the estimator is 


i 
(P11)st = S Wh(bn/ar) . (9) 
h=1 


A consistent estimator for the variance of (6 11)st is 


v(P11)st=_ Wh? v(bp/an) , (10) 
h=1 


where v(bh/ah) is defined in equation (8). 


The proportion of all sites occupied by a pair of owls was estimated separately for 
each subpopulation. Region-wide estimates were provided for each subpopulation. 
Estimates for each of the five physiographic provinces (table 1) were also provided 
for the nonreserved SOHA subpopulation. 


The stratified estimator (equation 5) for a population proportion was used to provide 
each estimate. The individual National Forests (table 1) were the strata used in 
estimating proportions (equation 5) and their standard errors (equation 6). Owing to 
small sample sizes in some National Forests, some forests were combined for 
computing standard errors. In all cases, National Forests that were combined were 
near each other geographically, so pooling was also supported by the geographic 
similarity of the units chosen to be pooled. Tables 2, 3, and 4 document the strata 
used to provide the estimates of the proportion of all sites occupied by a pair of owls 
and their associated standard errors. The size (number of potential sampling units) of 
each strata and the sample size in each strata are also provided. 


Table 2—Sizes of strata’ and samples in the nonreserved SOHA subpopulation 
used to estimate the proportion of sites occupied, by physiographic province 


Number of 
Physiographic National nonreserved Number of SOHAs 
province Forest SOHAs sampled° 
OLPE Olympic 25 25 
WACA Gifford Pinchot 35 14 
Mount Baker-Snoqualmie 44 31 
ORCA Deschutes 8 2 
Mount Hood 45 13 
Rogue River 16 3 
Umpqua 38 13 
Willamette 62 22 
Winema i 2 
CORA Siuslaw 25 24 
KLAM Rogue River 8 5 
Siskiyou 22 9 


? Strata are National Forests. 
Number is stratum size. 
° Number is sample size. 


Table 3—Sizes of strata? and samples in the reserved SOHA subpopulation 
used to estimate the proportion of sites occupied, by National Forest 


Number of Forest strata 
National Number of reserved SOHAS combined for 
Forest reserved SOHAs sampled° variance estimation 


Olympic 4 
Mount Baker-Snoqualmie 41 
Gifford Pinchot 7 
Mount Hood 6 
Willamette 16 
Siuslaw 15 
Rogue River® 4 
18 

3 

5 

2 

2 


< 


Siskiyou 
Umpqua 
Rogue River® 
Deschutes 
Winema 


— 
DONWOD$]$$H]ONDWMNDAW 
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* Strata are National Forests. 

Number is stratum size. 
° Number is sample size. 

National Forest strata followed by the same letter were combined for estimating the variance. 
° The Rogue River National Forest was divided into two strata. 


Table 4—Sizes of strata’ and samples in the nonreserved non-SOHA sub- 
population used to estimate the proportion of sites occupied, by National Forest 


Number of Number of Forest strata 
National sampling ynits sampling units combined for 
Forest on Forest sampled° variance estimation 
Olympic UUs 5 
Mount Baker-Snoqualmie 228 9 
Gifford Pinchot 259 9 
Mount Hood 157 5 
Willamette 263 10 
Siuslaw 122 5 
Rogue River 109 5 
Siskiyou 179 7 
Umpqua 184 8 
Deschutes 164 4 x 
Winema 26 1 x 


? Strata are National Forests. 
Number is stratum size. 
° Number is sample size. 
? National Forest strata followed by the same letter were combined for estimating the variance. 


Table 4 illustrates the combining of strata to provide an estimate of standard error. 
All 11 forests were considered as strata for estimating the proportion of occupied 
sites in the nonreserved non-SOHA subpopulation. The Winema National Forest 
strata had a sample size of one, so it was not possible to compute a standard error 
for that strata. The Deschutes and Winema National Forests therefore were com- 
bined and a standard error was computed for this combined strata to produce a 
combined strata with a sample size of five (4+1) and a strata size of 190 (164+26). 
This one combined strata was used in place of the two individual strata for computing 
s.e.(Pst) (equation 6). 


Occupancy by a pair of owls with young can be computed by using two different 
proportions that provide different pieces of information. One way is to compute the 
proportion of sites occupied by a reproductive pair that is based on all sites in the 
population (equation 5). A second method of computing the proportion of sites occu- 
pied by a reproductive pair is to base the proportion only on sites occupied by pairs. 
In this case, a random variable, the number of sites occupied by a pair, is the denom- 
inator of the computed proportion. An alternative estimator (equation 9), designed to 
estimate a proportion of a proportion, must be used in this case. The associated 
estimator of standard error is given in equation (10). Both methods were used so that 
the two different pieces of information provided by the two methods are available. 


In 1987, no data on occurrence of young were collected for the nonreserved non- 
SOHA subpopulation, so estimates of occupancy by pairs with young could not be 
provided for this subpopulation. For the reserved and nonreserved SOHA subpopula- 
tions, all estimates were provided only for these two land classifications combined 
but with the five original physiographic provinces collapsed into three groups of prov- 
inces. The three groups of provinces were (1) OLPE and WACA, (2) ORCA, and 


Table 5—Sizes of strata? and samples in the combined reserved and non- 
reserved SOHA population used to estimate the proportion of sites occupied by 
reproductive pairs, the proportion of occupied sites that contained reproductive 
pairs, and the average number of young produced over all sites, by physio- 
graphic province group 


Physio- Forest strata 
graphic Land Number of Number of combined for 
province classi- National areas in areas variance 
group fication Forest Forest sampled° — estimation 
OLPE Nonreserved Gifford Pinchot 35 14 
and Mount Baker-Snoqualmie 44 31 
WACA Olympic 25 25 
Reserved Gifford Pinchot U 6 
Mount Baker-Snoqualmie 41 18 
Olympic 4 3 
ORCA Nonreserved Deschutes 8 2 
Mount Hood 45 13 
Rogue River 16 3 
Umpqua 38 13 
Willamette 62 22 
Winema 7 2 
Reserved Mount Hood 6 3 
Rogue River 5 2 
Umpqua 3 3 
Willamette 16 6 
CORA Nonreserved Rogue River 8 5 
and Siskiyou 22 9 
KLAM Siuslaw 25 24 
Reserved Rogue River 4 1 4 
Siskiyou 18 6 Xx 
Siuslaw® 15 1 x 


? Strata are National Forests. 
Number is stratum size. 
° Number is sample size. 


? National Forest strata followed by the same letter were combined for estimating the variance. 
° The Siuslaw reserved strata was omitted from the population when the proportion of occupied sites 
containing a reproductive pair was estimated. The one sample unit in this strata was not occupied. 


(3) KLAM and CORA. Estimates were provided for these three groups of provinces, 
as well as for the entire Region, for the combined reserved and nonreserved sub- 
populations. Estimates for each of these situations were computed by using both 
methods of estimating occupancy by pairs with young as described above. The orig- 
inal strata were maintained for estimation; for example, typical strata are two from the 
Umpqua National Forest. The reserved SOHA lands and the nonreserved SOHA 
lands on the Umpqua National Forest are each separate strata. 


Table 5 indicates that the reserved and nonreserved SOHA subpopulations were 
combined to produce estimates of the proportion of all sites occupied by reproductive 


Estimating 
Reproduction 


10 


pairs for the physiographic province groups and for the Region. To produce esti- 
mates of the standard error, however, it was necessary to combine the reserved 
SOHAs of the Rogue River, Siskiyou, and Siuslaw National Forests in the CORA and 
KLAM physiographic province group. This was done because the Rogue River and 
Siuslaw National Forests each had only a single sample observed in the reserved 
land status. 


Estimates of the proportion of occupied sites that were occupied by a reproductive 
pair used the same strata as described above with one exception. The single excep- 
tion is the reserved SOHA lands in the Siuslaw National Forest that were omitted 
from the population. This occurred because the single sample observed was not 
occupied by a pair; thus, no information pertinent to this proportion for this National 
Forest was available. These strata are summarized in table 5. 


Estimates of the average number of young per site were provided for the three 
physiographic province groups (defined above) and the Region. The stratified estima- 
tors for the population mean (equation 1) and its standard error (equation 2) were 
used to obtain estimates for the combined reserved and nonreserved SOHA sub- 
populations. The strata were the same as previously described for estimating the 
proportion of sites occupied by a pair. 


In accordance with sampling protocol, sites were classified as having owl reproduc- 
tion based on observations of young. Because of the probability of predispersal mor- 
tality of young, we assumed that total number of young could be accurately counted 
only before July 1. We assumed that after July 1 it was possible to accurately deter- 
mine whether or not reproduction occurred at a site, but that late in the season, the 
observed number of young may not equal number of young produced. Therefore, for 
the seven sites where reproduction did occur but young were first observed after 
July 1, total number of young produced could not be accurately determined. For 
these sites, an overall average number of young, computed from all other sites 
where reproduction occurred, was used in calculations. This average was based on 
the 17 nest sites where reproduction was documented and number of young pro- 
duced was accurately determined before July 1. If the seven censored samples had 
been treated as unsampled units, an unacceptable bias would have been produced 
in the estimated number of young per site because only sites with observed reproduc- 
tion could be excluded from the sample. On the other hand, substituting an overall 
simple mean for these missing values may have produced a bias in the estimate of 
standard error. It is not clear, however, what the magnitude or direction of this bias 
was, because final estimates were in the context of stratified sampling. Further 
research is needed to decide the best approach for dealing with such missing values. 
One viable alternative is to use some method of statistical imputation (Sande 1982) 
to replace missing values with data. Sande (1982) provides a general discussion of 
methods of statistical imputation and problems with their implementation in practice. 


Estimates of the average number of young per site occupied by a pair were provided 
for the three physiographic province groups and the Region. The usual stratified 
estimators for the population mean (equation 1) and its standard error (equation 2) 
are appropriate if the sample from each strata is assumed to consist only of the 
occupied sites. However, the sampling fraction in stratum h, fh, is calculated based 


Sources of 
Potential Bias 


on the complete sample size and population size for the appropriate stratum, be- 
cause the population size for the occupied sites is unknown. The strata used for 
these estimates are the same as those described for the estimates of the proportion 
of occupied sites that were occupied by a reproductive pair. 


The field protocols used for making observations and measurements created sources 
of potential bias in some estimates. The protocols are briefly described in the “FY87 
Annual Report” (see footnote 1) with additional details in the “Draft Monitoring and 
Inventorying Handbook.” The main problem with the field protocol was that two 
types of measurement errors could have occurred and introduced bias into some 
estimates. These sources of potential bias are discussed below. 


One measurement error that could have occurred was misclassifying the occupancy 
of a site. A site occupied by a pair of owls may have been incorrectly classified as 
not occupied. Given the protocol procedures, the reasons for this misclassification 
are that the number of visits may have been inadequate or that some pairs may 
remain undetected regardless of the number of visits. In 1987, sites were visited a 
maximum of seven times. It may also have been possible to misclassify a site as 
occupied when, in fact, it was not occupied. This could occur if a male and female 
were observed on a site but were not actually paired. Specific precautions were 
taken in the protocol, however, to avoid this mistake. We therefore believe that the 
problem caused by the first kind of misclassification error is the greater source of 
potential bias, and it will be discussed in more detail. 


Classical sampling methods for estimating a proportion involve classifying each unit 
in the sample into one of two categories. The classification is assumed to be done 
without error, and few references discuss problems in estimation where errors in 
classification occur. Any such misclassification would result in biasing the proportion 
of sites that were occupied. 


The magnitude of bias caused by pairs undetected regardless of the number of visits 
is not known. There is a chance that if pairs were undetected, they might be quite 
different in traits from those detected with the established protocol. Perhaps undetect- 
ed owls are less likely to respond to calling or to offers of mice for food, or are differ- 
ent in any one of several other ways. This sampling methodology assumed that this 
is not the case. If this is the case, then the bias caused by misclassification would be 
greater than we suspect. This potential for bias warrants further investigation. 


Bias caused by misclassification because units are visited a maximum of only seven 
times may be reduced by two methods. One method is to increase the maximum 
number of revisits before a unit is classified as not occupied. This should reduce, at 
least theoretically, the number of misclassifications and hence any associated bias. 
Each revisit is quite expensive, however. The empirical frequency distributions of 


3 USDA Forest Service, Pacific Northwest Region. Draft 
monitoring and inventorying handbook. Rev. April 6, 1987. 
On file with: Fish and Wildlife Staff, Pacific Northwest Region, 
P.O. Box 3623, Portland, OR 97208-3623. 
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Figure 1—The observed frequency distributions of the number of visits required to establish 
occupancy. Distributions for the combined data from nonreserved and reserved SOHAs, as 
well as for the rionreserved non-SOHAs, are provided. 


number of revisits required to establish occupancy (fig. 1) indicates that only a small 
percentage of occupied sites were classified on the seventh visit. We feel this 
indicates additional revisits would probably not be a cost effective way to reduce 
misclassifications. 


A second approach to the misclassification problem associated with the number of 
visits is to statistically adjust the estimated proportion in an attempt to reduce the 
bias caused by misclassifications. An adjustment can be derived based on a prob- 
abilistic model of detecting the presence of an owl on a particular visit. Such an 
adjustment scheme was developed by Jim Baldwin.* The assumptions necessary for 
Baldwin’s statistical adjustment seemed very restrictive, however. The key assump- 
tions were that (1) sites are completely independent, (2) successive visits to the 
same site are independent events, and (3) the probability of detecting a pair of owls 
on any particular visit is a fixed constant. In particular, the probability of detection 
does not change over time or between visits. Serious questions were raised about 
the adequacy of the data in meeting the assumptions necessary for this adjustment 
procedure. Because of the problems of meeting the assumptions, the decision was 
made that no adjustment should be made at this time. If an improved adjustment 
procedure were developed or the assumptions necessary for Baldwin’s adjustment 
could be verified with data, such an adjustment might be used in future surveys. 


# Personal communication with Jim Baldwin, Pacific South- 
west Forest and Range Experiment Station, 1960 Addison 
St., Berkeley, CA 94704. 
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Spotted owls (Strix occidentalis) were monitored on 11 National Forests in the Pacific 
Northwest Region of the USDA Forest Service between March and August of 1987. 
The basic intent of monitoring was to provide estimates of occupancy and reproduc- 
tion rates for pairs of spotted owls. This paper documents the technical details of the 


statistical estimation procedures used to produce all the estimates resulting from that 
effort. All relevant formulae are presented as well as decisions made about collap- 
sing strata to estimate variances. Some complicating aspects of the monitoring effort 
are also discussed. 
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